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We explore the impact of the short range interaction on the scattering of ground state polar 
molecules, and study the transition from a weak to strong dipolar scattering over an experimentally 
reasonable range of energies and electric field values. In the strong dipolar limit, the scattering scales 
with respect to a dimensionless quantity defined by mass, induced dipole moment, and collision 
energy. The scaling has implications for all quantum mechanical dipolar scattering, and therefore 
this universal dipolar scaling provides estimates of scattering cross sections for any dipolar system. 
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Collisional control in ultracold atomic physics has led 
to the study of many remarkable systems such as ultra- 
cold molecules [ij and the BEC-BCS crossover Q. Re- 
cently, collisional control of chromium has been achieved 
, resulting in a clear demonstration of magnetic dipole- 
dipole interaction Q . Other experiments are rapidly pro- 
gressing towards the production of ground state polar 
molecules (GPMs), which have large permanent electric 
dipole moments. These interact via the dipole-dipole in- 
teraction: 



Ai • A2 - 3(R ■ Ai)(R ■ A2) 
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where A is the electric dipole moment of a molecule and 
R is the intermolecular separation. Because this interac- 
tion is long range and anisotropic, its incorporation into 
many body systems has led to exciting predictions such 
as dipolar crystals in restricted dimensions [H] and ro- 
tons Q . There are also many intriguing applications 
of GPMs, such as quantum computing [1, Q. 

Many different techniques are being used to obtain 
GPMs [13] , most notably is photo-association (PA) . The 
PA experiments are nearing production of ultracold polar 
molecules in their absolute ground states for a ran ge o f 
different heteronuclear alkali- metal systems ^,12, 13| . 

molecules are relatively simple; only the rotational 
structure needs to be accounted for to study the colli- 
sions accurately. 

To set the scene, first consider ultracold atomic sys- 
tems, where collisions are parametrized by the scatter- 
ing length, as- At a Feshbach resonance the scatter- 
ing threshold becomes degenerate with a molecular state, 
leading to a divergence of as . Thus the magnetic Fesh- 
bach resonance allows tuning of the interaction, but this 
is solely a short range affair. Furthermore the only impor- 
tant scattering is s-wave, unless there is another resonant 



partial wave, e.g., see Ref. [ij]. This scenario is in direct 



contrast to dipolar scattering as will be shown. 

For weakly dipolar systems, such as chromium, as 
plays a significant role in determining the dynamics of 
the system. Chromium experiments are now exploiting 



Feshbach resonances to tune as near zero; so dipolar in- 
teractions are dominant, and in some cases leading to 
dipolar collapse In contrast, for strong dipolar scat- 
tering the short range interactions play a minor role, if 
not irrelevant. 

In anticipation of these molecular systems, we have 
studied the scattering of GPMs over a wide range of 
electric fields and collision energies to determine the in- 
fluence of the short range interaction. This work shows 
the scattering of dipoles is only weakly influenced by the 
short range interaction. We also see the emergence of a 
universal scaling of dipolar collisions which has implica- 
tions for all quantum mechanical scattering dipoles. 

To understand Eq. ([T]) more clearly, consider the in- 
teraction in its asymptotic form in the lowest threshold 
containing two GPMs. This is achieved when R is large, 
and so couplings to higher thresholds are negligible. This 
distance is quite large, typically greater than lOOoo for 
GPMs, where oq is the Bohr radius. The long range in- 
teraction is proportional to the induced dipole moment, 
and this requires a non zero electric field. For simplicity 
of notation we use d to denote the induced dipole moment 
created by a field along the z axis; it is the expectation 
value of the dipole moment of the field dressed ground 
state. With these assumptions, the dipole-dipole interac- 
tion has the form: (l — 3 cos^(6')) , where 9 is the angle 
between R and z. This illustrates the anisotropy clearly 
and shows that the strength of the interaction relies di- 
rectly on the induced dipole moment. The anisotropy 
and long range nature of the interaction induce signifi- 
cant scattering contributions from non-zero partial waves 
Furthermore, the long range nature of the in- 
teraction alters threshold behavior. The Born approxi- 
mation predicts that all partial wave cross sections are 
constant as the collision energy goes to zero 17 1. 

The characteristic length and energy of a dipolar sys- 
tem are defined in terms of mass m (or twice the reduced 
mass) and the induced dipole moment d, and they are 
D — md"^ and Ed — d'^/D'^ = 1/irt'd'^, respectively. 
With the scattering energy, we can form a dimensionless 
quantity ^ — E/Ed — m^d'^E which parametrizes the 
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scattering. At a fixed field it is the dimensionless en- 
ergy of the system. A characteristic electric field for a 
molecule, Sq, is determined hy B/ji where B is the 
rotational constant of the GPMs. 

To study the GPM systems we determine the total 
cross section cr, which is 
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where ri/^-* is the T matrix which details the collisions 

''J 

leading to a transfer from channel i to j for a system with 
azimuthal symmetry M about the field axis 18| . k is the 
wave number, \/raE. The factor of two in Eq. ^ is there 
for only initially identical scatters. In this paper we as- 
sume identical bosons unless stated otherwise. We per- 
form scattering calculations for several polar molecules 
to obtain T and a. These are extremely large computa- 
tional tasks due to the large number of partial waves and 
total Ms required to converge the calculation because of 
the anisotropy and long ran ge n ature of Eq. ([T]). The 
details are presented in Ref. [l6|. Here we have added a 
Lennard-Jones potential and vary the inner wall to alter 
the zero field scattering length. The minimum of this 
short range potential, Rmim is typically 10 ao and much 
deeper than V^^(-Rmm), thus changing the character of 
the scattering potential. In zero field {d = 0) the GPMs 
scatter similarly to atoms and are parametrized by as- 

An estimate of the total cross section is achieved with 
a semi-classical approach. This approach offers scaling 
of (7 on the physical parameters of the system, such as 
d, m and E jiol . |20| . Another important cross section is 
the quantum unitarity limit, which provides a maximum 
value for any single partial cross section. This occurs 
when the T matrix takes on its maximum value of 4. 
These cross sections are: 
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where csc — 1-7x10^^^ and cq = 4.85 x 10^^^ are chosen 
so that the units of /i, m, E, and a are [D], [a.m.u.], [K], 
and [cm^], respectively. 

To begin the analysis we illustrate the influence of the 
short range interaction as the electric field is varied. We 
have calculated a of RbCs 2^ as a function of energy for 
many different electric fields and with different values of 
Qs- Figure[T](a) shows the total RbCs cross section versus 
energy at electric fields of 0.5 (black), 1 (red), 2 (blue), 3 
(green), and 5 (brown) Eq . For each electric field there 
are 7 different values of Os shown ranging from -650 to 
500 ao. The difference between these calculations at a 
given electric field is only the inner wall coefficient, and 




10 10 

Energy (K) 




10 10 
Energy (K) 



FIG. 1: (Color Online) (a) The variation of the cross section 
for many different electric fields with different short range po- 
tentials. Different fields correspond to sets of colored curves; 
in ascending order the fields are 0.5 (black), 1 (red), 2 (blue), 
and 5 (brown) So- (b) The average total cross section, a, 
(black) for various field values (in ascending order 0.5, 0.75, 
1, 1.5, 2, 3, 4, and 5 £o) and the variations, Aa, about a (blue 
with circles). In both plots the red dashed line is ctq and the 
green dashed-dot line is proportional to E^^^^ as suggested 
by asc- 



the calculations are identical for R > 



The varia- 



tion between calculations is due solely to the difference 
in short range potential, and therefore these calculations 
directly access the influence of the short range potential 
on the dipolar scattering. The essential result of this flg- 
ure is that at high electric field and high energy or large 
^ (brown) the scattering is insensitive to the short range 
interaction. This is in contrast to low electric field and 
low energy or small ^ (black), where there is great varia- 
tion in the scattering due to different phases acquired at 
short range. 

To study this figure in detail, first we look at the weak 
field results (black), small ^. The GPMs are able to ac- 
cess the short range interaction and therefore scattering 
is sensitive to this process. The cross section can be res- 
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onantly large because erg > crsc [ifil , and these resonant 
variations can dominate the scattering. As the electric 
field is increased to 1 (red) and 2 (blue) £q, the cross 
section becomes larger and the variation in the scatter- 
ing cross section becomes relatively small; this is most 
evident at high energy. Finally, at a large electric field, 
5 Eq (brown), the cross section is very large and there is 
only slight variation. The dipole-dipole interaction has 
induced large numbers of non-zero partial waves to the 
scattering; we note a >> uq. These non-zero partial 
waves are insensitive to the short range interaction. This 
fact constrains the resonant control of the scattering as 
has been seen in ultracold atoms with magnetic Fesh- 
bach resonances. It is also important to note that the 
cross sections do not go to zero even for small fields, for 
similar reasons. 

To extent the analysis, we have averaged a at each 
field (ct) shown in [1] (b) . In ascending order the fields 
for a are 0.5, 0.75, 1, 2, 3, 4, and 5 £o (black). We have 
also obtained the variation of the total cross section Act, 
Acr^ — Y^llii^a — o'Y ■ The variations in a are shown 
in Fig. [I](b) as blue full circles. The largest variations oc- 
cur from resonant scattering of a partial wave, and have 
a maximum contribution of aq . So it is reasonable to ex- 
pect Act ~ aq. This leads to an interesting comparison, 
in the large ^ limit, which is 
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FIG. 2: (Color Online) The transition of dipolar scattering to 
a universal behavior is shown by plotting T vs ^ for many po- 
lar molecules. The molecules are *^Rb^^K (black x), fermionic 

The 
The 

blue line with triangles is proportional to ^. The inset is the 
experimental cross section for Rydberg atoms from Ref [20| . 
The blue line with open circles is the thermally averaged Eq. 
0- 



^^^Rb^'^'K (black +), NaCs (brown □), and RbCs (red o) 
unitarity limit is 4, and is denoted by the dashed line 



Act ctq 1 
^ — — cx -7=. 



(5) 



This offers a more explicit statement of why for large ^ 
there is minimal influence of the short range interaction. 

Once the short range interaction is negligible compared 
to the long range interaction (large ^), it is instructive 
to rescale the multi-channel radial Schrodinger equation 
using the length scale D. Including only the long range 
influences of kinetic energy and the dipole-dipole inter- 
action, one finds that the only parameter left in the set 
of equations is ^. Performing the rescaling we find 
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where y = r/D and ipf^ is a multi-channel radial wave- 
function. The coupling between partial waves, C^^f, is 
well known [2^. Since Eq. ^ only depends on ^, it 
suggests universal scaling of dipolar scattering. To il- 
lustrate this behavior we have compiled scattering data 
from many different molecular systems. In Fig. [2] we 
have plotted T as a function of ^ for ^^Rb'^^K (black 
x), fermionic ^^Rb^°K (black -I-), NaCs (brown squares), 
and RbCs (red diamonds), which is the data in Fig. [T] 
(a) [i^l . A line proportional to ^ is shown as a blue line 
with triangles. This line represents a cross section which 
is constant as ii^ ^ at a fixed field. The figure shows 
the transition in the scattering from highly variable at 



low ^ to uniform at large ^. This transition of T sig- 
nifies the onset of universal dipolar behavior. This will 
occur when the dipolar interaction is dominant and the 
scattering will be insensitive to the short range interac- 
tion. For this reason different molecules, even bosons and 
fermions, have the same scattering behavior. 

In Eq. (I4|), erg is defined by the T matrix taking on 
its maximum value of 4. For the same reason as Ac 
aq, we find AT ^ 4 for all scattering. For small ^, 
the scattering can access the short range and therefore 
resonant scattering is significant and so are the details of 
the short range. This is seen by the great variation in the 
scattering data for ^ < 200. The transition to universal 
scattering behavior is seen as ^ is increased above 200. 
The possible values of T initially span a wide range, but 
this span greatly decreases at high This is due to 
large contributions from many non-zero partial waves; so 
typical values of T are much greater than the variations 
of any single term. 

The universal dipolar behavior is clearly seen as the 
scattering of the dipoles becomes uniform at high ^. Fit- 
ting T at large ^, we find 

r = 0.266Vi- (7) 

This is shown in Fig. [2] as a blue line with full circles. 
This result was obtained by fitting all RbCs and NaCs 
data from ^ — 10^ to 10^, i.e., 458 scattering calculations. 
This simple equation offers an estimate of T and a for 
all quantum mechanical scattering dipoles. 
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For large ^, dipolar systems obey a universal scaling, 
where all scattering dipoles will behave similarly irre- 
spective of the details of the short range. This implies 
the dipoles can be bosons, fermions, identical, or distin- 
guishable, and the theory will apply. A striking example 
of this theory being applied is an experimental measure- 
ment of the cross section for resonant collisions of Ryd- 
berg atoms [20|. In this experiment two identical Ryd- 
berg atoms in the ns state, where n (s) is the principal 
quantum number (orbital angular momentum), are reso- 
nantly scattered into a degenerate threshold to which it is 
coupled via the dipole-dipole interaction. For a particular 
electric field, the ns + ns threshold becomes degenerate 
with the np + (n — l)p threshold. This system has huge 
dipole moments d cx n^, e.g., consider n = 22, the dipole 
moment is about 100 D\ Numerically converging this 
calculation would be impossible with the present compu- 
tational techniques. But using the scaling presented here 
we can obtain an accurate estimate of the total cross sec- 
tion. 

The inset in Fig. [2] contains the experimental cross sec- 
tion (squares) for sodium Rydberg atoms [20]. The blue 
line with open circles is Eq. ([7]) with an additional factor 
of y^Tr/W to account for thermal averaging of collisions 
in a beam [2j] . To calculate the values of ^ for the exper- 
imental data we use d = 0.6n*^, where n* is the effective 
quantum number, and an average coUisional velocity of 
w = 1.6 • 10"'' a.u. (T = fv^ ~ 0.17i^). The agreement 
of the slope is not surprising [l^, [2^ , but agreement of 
the magnitude is quite good. There has never been a 
means to determine the amplitude with accuracy. The 
coefficient in Eq. ([7]) allows us to accurately predict the 
cross sections for all scattering dipoles. The limit of this 
theory is when other physics emerges and alters the form 
of Eq. ([6|). For example, in many-body systems there will 
be screening and other molecules will become important 
before the two scattered molecules leave the scattering 
volume. 

We have studied the scattering of ultracold ground 
state polar molecules for experimentally accessible en- 
ergies and fields. We have illustrated the limited in- 
fluence of the short-range interaction in the presence of 
dipolar interactions. We have also rescaled many dif- 
ferent scattering calculations, finding dipolar collisions 
are parametrized by ^ = nr'd'^E. For our calculations ^ 
ranges up to 10^. This is well into the region of universal 
dipolar scaling, and therefore we were able to determine 
Eq. ([7]) with accuracy. This equation predicts collision 
cross sections for all scattering dipoles. Consequently 
we are able to estimate cross sections for quantum me- 
chanical dipolar systems which are far beyond our cur- 
rent computational capabilities. Future directions of this 
work will be to study the effects of nearly degenerate 



thresholds on the scattering, such as those presented by 
hyperfine structure. 
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